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Electroweak bosons produced in conjunction with jets in high-energy collider experiments is one 
of the principle final-state channels that can be used to test the accuracy of perturbative Quantum 
Chromodynamics calculations and to assess the potential to uncover new physics through compari- 
son between data and theory. In this paper we present results for the ^'^/7*-|-jet production cross 
sections at the LHC at leading and next-to-leading orders. In proton-proton reactions we elucidate 
up to 0{Gf(xI) the constraints that jet tagging via the Z°/7* decay dileptons provides on the 
momentum distribution of jets. In nucleus-nucleus reactions we demonstrate that tagged jets can 
probe important aspects of the dynamics of quark and gluon propagation in hot and dense nuclear 
matter and characterize the properties of the medium-induced parton showers in ways not possible 
with more inclusive measurements. Finally, we present specific predictions for the anticipated sup- 
pression of the Z° /^* +jet production cross section in the quark-gluon plasma that is expected to be 
created in central lead-lead collisions at the LHC relative to the naive superposition of independent 
nucleon-nucleon scatterings. 

PACS numbers: 12.38.Bx, 13.87.-a, 12.38. Mh 
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I. INTRODUCTION 

Hadronic jets [l| produced in today's high-energy col- 
lider experiments have long been regarded as a premier 
tool to test the fundamentals of perturbative {Q^ ^ 
^qcd) Quantum Chromodynamics (QCD) Q. The 
start-up of the Large Hadron Collider (LHC) has stim- 
ulated new theoretical and experimental interest in jet 
observables @. At the LHC, final states involving 
high jet multiplicities, jets-|-long-lived heavy leptons and 
jets-l- missing energy are among the most-promising chan- 
nels for discovery of physics beyond the Standard Model. 
It is, therefore, critical to understand theoretically as ac- 
curately as possible the QCD background to signatures 
of new physics. 

One of the most studied QCD final states at collider 
energies is an electroweak boson accompanied by jets [3- 
0]. Of these processes, the Z"/7*+jet production is rel- 
atively easy to measure via the Z°/7* + l~ de- 
cay channel. Comparison between theory and experi- 
ment has so far only been carried out at the Tevatron 
-y/s = 1.96 TeV [1] without restrictions on the momen- 
tum of the vector boson. In this paper we elucidate the 
constraints on the production and momentum distribu- 
tion of jets associated with Z^/j* of fixed pr to low- 
est [OiGpas)] and ncxt-to-leading [0{GFa1)] orders 0- 
The cross sections that we find are also the much needed 
baseline for the investigation of many-body QCD effects 
in heavy ion reactions at the LHC. 



A significant part of the LHC research program in- 
volves reactions with ultrarelativistic nuclei - specifically, 
Pb-|-Pb (lead-lead) collisions up to a center-of-mass en- 
ergy per nucleon pair ^Jsnn = 5.5 TeV. These heavy 
ion runs will consolidate the evidence for the creation 
of a deconfined state-of-matter at extremely high tem- 
perature, the quark-gluon plasma (QGP), and probe its 
properties in a new energy regime [9| . One of the prin- 
cipal advantages of the LHC over the Relativistic Heavy 
Ion Collider (RHIC) is that the projected 28-fold increase 
in energy will open final-state channels for jet tomogra- 
phy of the QGP that are currently inaccessible to ex- 
periment or limited by statistics [10[. One such channel 
is Z°/7*-|-jct, which has been previously discussed in the 
context of the energy or momentum constraints that tag- 
ging of the decaydileptons can provide on the away-side 
parton shower [ll| - [l3j . The short production time, 
1/to^, and subsequent decay, Tdecay ^ l/P, im- 
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ply that these processes occur before or during the for- 
mation of the QGP. While the produced dileptons have 
to traverse a region of dense nuclear matter, they reach 
the detectors unscathed by the strong interactions in the 
deconfined QCD medium. 

The utility of jet tagging, however, extends well be- 
yond the simple determination of the energy of the par- 
ent parton. It allows one to probe aspects of the QCD 
many-body dynamics in nucleus-nucleus (A-l-A) colli- 
sions that are inaccessible via leading particles or even 
leading particle correlation measurements [l^ |. One of 
the most striking results from the heavy-ion program at 
RHIC [l^, [l^ is the attenuation of the flux of energetic 
hadrons in the plasma, known as jet quenching. Much 
theoretical effort has been invested in understanding the 
mechanism of this suppression a nd, in p articular, of par- 
tonic energy loss in the QGP |17H2(]| . Unfortunately, 
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FIG. 1: (Color online) Experimental results from the Fermi Lab Tevatron Collider [8| for the cross section of jets associated 
with Z'^ /■y* — >■ /i+ + M~ in p + p collisions at ^/s = 1.96 TeV are compared to tree level pQCD. The left panel shows this cross 
section versus pj, q^jj, whereas the right panel shows it versus J/(jct)- The (red) lines indicate our LO result. The (yellow) bands 
indicate the variation from the choice of scale. The experimental cuts are discussed in the text. 



leading particle quenching alone [21[ is not sufficient to 
discriminate between partonic energy loss formalisms or 
to extract quantitatively the stopping power of the QGP 
for color-charged particles [l^l- In contrast, jet obscrv- 
ables arc much more closely related to the underlying per- 
turbative QCD theory and to the characteristics of the 
vacuum and medium-modified parton showers. Hence, 
they are also much more discriminating with respect to 
theoretical approximations and model assumptions. 

In this paper, we extend the formalism developed to 
evaluate inclusive jet cross sections and shapes in A-t-A 
reactions at lowest and next-to-leading orders [lol . fl^ to 
tagged jets. We present the results of a systematic study 
of the production and subsequent suppression of the 
Z'^/7*+jct final state in central Pb-|-Pb collisions at the 
LHC at y/sWN = 4 TeV to OiGpal). Wc demonstrate 
that tagged jets can provide insight into the multiplicity 
and distribution of gluons induced by parton propagation 
in the QGP. Much of the discrepancy in the determina- 
tion of the QGP properties can be traced to differences 
in the characteristics of these distributions [1^ [l^l- In 
our complete numerical examples for the modification of 
the /'y* tagged jet cross section we use the Gyulassy- 
Lcvai-Vitev (GLV) approach to the non-Abelian energy 
loss of quarks and gluons propagating in dense nuclear 
matter [25l |. 

Our work is organized as follows: in section II we 
present results for the cross section of the /j* + jet 
channel at the Tevatron and at the LHC. Wc elucidate 
the crucial differences between the lowest-order and ncxt- 
to-leading order calculations for the py-tagged jet cross 
section where the transverse momentum of the boson is 
constrained in a narrow interval via its decay products 
In section HI we discuss one of the principle dif- 
ferences between proton-proton and nucleus-nucleus col- 
lisions: the radiative corrections from final-state interac- 
tions of Z°/7*-taggcd quark and gluon jets in the QGP. 
Wc give details for the calculation of the medium in- 



duced bremsstrahlung at LHC energies. Our results for 
the suppression of the tagged jet cross sections in cen- 
tral lead-lead (Pb+Pb) collisions are given in section IV. 
We demonstrate how tagged cross sections can probe es- 
sential aspects of the physics of parton energy loss in 
strongly-interacting matter in ways not possible with in- 
clusive jets or leading particles. Our conclusions are given 
in section V. Appendix A contains selected steps in the 
calculation of the lowest order tagged cross sections. In 
appendix B we briefiy discuss the Dalitz decays of the 
Z'^/j* to dileptons. Finally, the sensitivity of the ex- 
perimentally observed cross section suppression to the 
fraction of quark- and gluon-initiatcd jets is discussed in 
Appendix C. 



II. Z°/7*-TAGGED JETS IN HADRONIC 
COLLISIONS 

We begin by discussing the cross section for Z°/j*- 
tagged jet production in hadron-hadron collisions. It is 
instructive to first consider the leading order (LO) re- 
sult, from which one can understand the underlying pro- 
duction processes and appreciate why the Z*^ boson was 
originally considered as a suitable tag for the initial as- 
sociated jet energy Details of this calculation are 
given in appendix El 

In the coUinear factorization approach, the Z*'/7*+jet 
cross section reads: 



da 



dy{z) dy(jet) <PpT{z)<PpT{ict) 



|M| 



/r, ^2 ^ - - C2 ^ (PT(Z) - PTQo 



(1) 



In Eq. ([T} S* = {pi + is the squared ccntcr-of-mass 
energy, f(xi,fi) arc the parton distribution functions, 
xi = I , xi — P2 / P2 the initial-state parton 
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FIG. 2: (Color online) Same as in Figure [5] but with a comparison of the Tevatron data to the MCFM LO result (blue) dotted 
line and the MCFM NLO result (green) solid line. Note that experimental error bars are statistical only. 



momentum fractions that are fully determined at tree 
level, and |Mp are the relevant squared matrix elements. 
The constraint on the transverse momentum of the jet is 
exact only at this order and only at the partonic level. 

As mentioned briefly above, in order to make connec- 
tion to experiment we will focus on the leptonic decay 
products of the Z° boson. For this reason, one must also 
include the contribution of virtual photon-tagged jets in 
the invariant mass range around the Z*^ peak. In our ap- 
plication to p+p and A-j-A collisions at the LHC we con- 
sider dilcpton pairs in the invariant mass rangcroziSr^, 
where mz = 91.2 GeV and F;, = 2.5 GeV U^. In this 
invariant mass range, the contribution from the dom- 
inates that of the virtual photon by roughly two orders 
of magnitude. To evaluate the tagged cross section in 
the dilepton decay channel, we use a Monte Carlo simu- 
lation to generate the isotropic fi'^ + distribution in 
the Z°/7* rest frame and boost this distribution back to 
the laboratory frame. Branching ratios arc taken from 
Ref. [2^ . Our approach allows us to precisely match the 
kinematic detector acceptance cuts. Further details of 
the Dalitz decay implementation are given in appendixlB] 

A comparison of our LO result for jets associated with 
Z°/7* n+ + in p4-p collisions at = 1.96 TeV 
with experimental measurements at the Fermi Lab Teva- 
tron Collider Q are shown in Figure [T] The experi- 
mentally measured muon pairs are in an invariant mass 
range of 65 GeV - 115 GeV. Additionally, the muons 
were required to have > 15 GeV and rapidity \y\ < 
1.7. Jets were experimentally reconstructed using a mid- 
point cone algorithm with Lorcntz-invariant cone size 
R = \J (A0)^ -I- (A?;)^ = 0.5. Furthermore, jets were 
required to have pt > 20 GeV and \y\ < 2.8. The 
tree level calculation under predicts the magnitude of the 
cross section by about 30% but describes well its shape 
both versus transverse momentum and rapidity. The 
(yellow) band illustrates the sensitivity of the Z^/'y*+iet 
cross section to the standard fi = [rriT /2, 2mT) variation 



of the renormalization and factorization scales, where 
77iT = \/ m\ -I- ]}\.. It can be taken only as a rough indi- 
cator of the size of one-loop corrections. 

At ncxt-to-leading order in as the Z°/7*-f jet cross sec- 
tion can be written schematically as: 

da 



dy(z) dyQct) <Pvt {z)(Pvt (jet) 

l_ 

dy(z) rfy(jct) <Pl*T {z)d'^PT (jet) 2! 



dy[z) dy{jot) dy3d'^pT{z)d'^PT(jct)d'^PTi 



(2) 



The two-particle final state represents the tree-level con- 
tribution and NLO virtual corrections. In the three- 
particle final state one of the jets is integrated out. For 
the NLO cross section evaluation we make use of the 
publicly available Monte Carlo for FeMtobarn processes 
(MCFM) 01, (available online at|http://mcfm.fnal.gov/ 1. 



MCFM provides one-loop results for many QCD pro- 
cesses of interest to the hadron collider physics commu- 
nity. Comparison between the MCFM predictions to LO 
and to NLO and the DO Tevatron II data is given in Fig- 
ure [21 At tree level the calculated cross sections are iden- 
tical to our results, up to differences between the CTEQ 
and MSTW parton distribution functions [H [H]. At 
NLO, the calculation is noticeably closer to the experi- 
mental data and within the systematic uncertainty of the 
measurement Q (not shown in Figures [T] and [5]) . 

The main advantage of the ncxt-to-leading order 
Z''/7*-|-jet calculation is the ability to precisely predict 
the transverse momentum distribution of jets associated 
with a dimuon tag in a narrow pt interval. Beyond tree 
level, the momentum constraint that the Z° boson mea- 
surement provides is compromised by parton splitting 
and Z-strahlung processes. We demonstrate this in Fig- 
ure [31 which shows the single differential cross section 
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FIG. 3: (Color online) Transverse momentum distributions of a jet associated with Z''/7*(— s> fi'^ + A* ) ^^S to 

and OlGpcts). Left and right panels show results for two transverse momentum cuts 92.5 GeV < pr < 112.5 GeV and 

7 GeV < pt < 13 GeV on the tagging particle, respectively. 



for jets tagged with /j* ji^ fJ,^ in p+p collisions at 
y/sNN = 4 TeV. Our choice of center-of-mass energy is 
motivated by the anticipated capabilities of the heavy- 
ion program at the LHC in the fall of 2010 We 
implement acceptance cuts of \y\ < 2.5 for both jets and 
final-state muons, and, as mentioned above, constrain the 
invariant mass of the muon pair to the interval m^iSFz, 
where mz = 91.2 GeV and = 2.5 GeV. This is the 
kinematic acceptance range in which we will evaluate all 
results that follow. 

For the cross section shown in Figure [3l the tagging 
Z77* is required to have 92.5 GeV < pr < 112.5 GeV 
(left panel) or 7 GeV < pt < 13 GeV (right panel). The 
LO result restricts the pr of the jet to lie exactly within 
this interval, consistent with Eq. ([1]). As seen in Figure |3l 
at NLO the deviations from this naive relation are very 
significant. We have included results for three different 
values of the jet cone radius, R = 0.2, 0.4, 0.8. The vari- 
ation of the cross section with R around Pt (jot) ^ Pt (Z) 
arises from the interplay between the amount of energy 
that is contained in the jet and the number of recon- 
structed jets (one or two). For PT(jet) ^ Pt(z) or 
PT(jet) ^ Pt(z) the two final-state partons are well- 
separated and identified as different jets. The falloff of 
the differential cross section relative to its peak value at 
Pt (jet) = Pt (z) is then controlled by the QCD splitting 
kernel (the part related to the large lightcone parton mo- 
mentum) and there is no dependence on the cone radius. 



In order to quantify the inability of the Z°/7* tag 
to constrain the momentum of the jet we calculate the 
mean p^ = {pt {jot}) ^J^d standard deviation Apt^q^j) = 

^{p't (jet)) ~ (PT(jot))^ for each of the curves in Figure [S] 
The results are presented in Table H] The standard de- 
viation for the LO curves is not strictly zero because of 
the finite pr width of the tagging Z'^/j* bin. The NLO 
curves exhibit a similar (pT(jot)) f-s the LO result, with 
(PT(jet)) increasing as the cone radius increases. How- 
ever, in going from LO to NLO there is a significant 
jump in Aprp Qet) ■ The width of the jet momentum distri- 
bution quadruples for the more energetic tag. The very 
large values of I^Pt (]ct) I {Pt {ict)) ^ 25% at NLO cre- 
ate serious complications for experimentally tagging the 
initial associated jet energy in both p-l-p and A+A colli- 
sions. While additional cuts can be considered, such as 
the requirement for a single jet within the experimental 
acceptance that is exactly opposite the tagging particle 
in azimuth, these will reduce the already small projected 
multiplicity for this final state in heavy ion reactions. 



III. MEDIUM-INDUCED PARTON SHOWERS 
FROM FINAL-STATE INTERACTIONS IN A-|-A 
REACTIONS 



A. Radiative energy loss of fast partons in the 
QGP 



The principle difference between jet physics in proton- 
TABLE I: Mean pr and A(pt) for Z°/7*-tagged jets at the proton collisions and jet-physics in nucleus-nucleus colli- 

LHC. sions is the contribution of final-state inelastic quark and 

gluon interactions in the QGP to the formation of parton 
showers [l^, [3 HO, HH • In this work we use a theoreti- 
cal approach developed to address the problem of parton 
energy loss at RHIC and at the LHC |l9l.[25t. 

To evaluate the medium-induced bremsstrahlung, it is 
important to keep track of the evolution of the gluon 



Pt(z) [GeV] 


LO 


R = 0.2 R 


= 0.4 R = 0.8 


7-13 (pt (jot)) [GeV] 


9.3 


8.7 


9.4 9.8 


Apt (jet) [GeV] 


2.2 


4.5 


4.7 4.6 


92.5-112.5 (pt (jot)) [GeV] 


100.8 


93.9 


96.6 100.1 


Apt (jot) [GeV] 


7.0 


25.2 


24.9 24.2 



5 



transverse momentum k in a plane perpendicular to the 
direction of jet propagation. Such a k may arise from 
a single hard scattering or from multiple soft scatter- 
ings. The acceleration of the color charges in the 2D 
transverse plane generates color currents whose detailed 
interference pattern determines the strength of the non- 
Abelian Landau- Pomcranchuk-Migdal effect [l^ [T7l - l20| . 
Let us denote by: 



B 



H 

(il---im)(jU2---in) 



k 

k2 ' 



(k - q,:i qi„)2 

H — C,;, , 



= c 



(3) 

(4) 

(5) 
(6) 



the Hard, Cascade, and Bertsch-Gunion propagators in 
the transverse momentum space [11. In Eqs. (0)-® q, 
are the momentum transfers from the medium. Another 



important quantity, which enters the bremsstrahlung 
spectrum, is the formation time of the gluon, t/, at the 
radiation vertex. When compared to the separation be- 
tween the scattering centers Azj ~ Zj — Zj-i, which can 
fluctuate from zero up to the size of the medium L, it 
determines the degree of coherence present in the multi- 
ple scattering process. We introduce the following nota- 
tion [ic|: 



k^ 

k+ ' 

k+ 
(k - q^i - 



(7) 
(8) 
(9) 



where is the gluon's large lightcone momentum. 

For final-state (FS) interactions, the double differential 
distribution of medium- induced gluons reads [25| : 



_ dN3{FS) ^ CrOs 
dk+(P\i 7r2 



E 



n / ^'^^ 

.7 = 1 \ 



1 daei{zj) 



m—l 



cos ^ 

^{k,--- ,n) Azk -cos y^tj(fc^...^„)Azfc 



.(10) 



, k=2 



In Eq. (Uni) = and B 



(n+l,n) 



B„ is under- 



stood. In the case of final-state interactions, zo ~ 
is the point of the initial hard scattering and z^ ~ L 
is the extent of the medium. The path ordering of the 
interaction points, zl > Zj+i > Zj > zq, leads to the 
constraint X^iLi ^-^^ — ^l- One implementation of this 
condition would be Azi £ [0, zl — J2]=^i ] and it is 
implicit in Eq. ([10]). We also note that Ci^ = Ci^ = 4/3 
for quark jets and Cr — Ca = 3 for gluon jets. The 
transverse momentum transfers from the medium are 
averaged over the normalized scattering cross section 
[l/aei{zj)]daei{zj)/d^qj. 



B. Numerical methods and QGP properties 

Results relevant to the LHC phenomenology are cal- 
culated using full numerical evaluation of the medium- 
induced contribution to the parton showers. These affect 
the observed jet shapes and the in-medium jet cross sec- 
tions. Energetic inclusive jet production and tagged jet 
production are rare processes that follow binary collision 
scaling ~ d^TVbin/rf^x. In contrast, the medium is dis- 
tributed according to the number of participants density 



^ o?^7Vpart/d^x. Soft particles that carry practically all 
of the energy deposited in the fireball of a heavy ion col- 
lision cannot deviate a much from such scaling. We take 
into account the longitudinal Bjorken expansion of the 
QGP and the medium density as a function of proper 
time reads: 



p(x,ro)— . 
r 



(11) 



In Eq. (|lip the initial density at time tq can be related 
to the experimentally measured charged particle rapidity 
density and the participant density in a plane perpendic- 
ular to the collision axis: 



p(x,To) 



1 d^{dN3/dy) 



To 



d2x 



1 3 
^2 



drj 
dy 



d^idN'^^/dr]) 



d2x 



d^N, 



part 



(12) 



As an illustration, we show in the left panel of Figure |3] 
the distribution of jets in central (& = 3 fm) Pb-|-Pb 
collisions at the LHC. This binary collision density is 
evaluated using an optical Glauber model with an in- 
elastic scattering cross section cri„ = 65 mb. The den- 
sity of QGP partons/soft particles in the plane transverse 
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FIG. 4: (Color online) The left panel shows the binary collision density in central Pb+Pb collisions at the LHC from an optical 
Glauber model with (Ji„ — 65 mb. The right panel shows the transverse density of soft partons/particles, which is proportional 
to the participant density. 



to the collision axis is shown in the right panel of Fig- 
ure m We note that for the constituents of the medium 
we use parton/hadron duality, sec Eq. ([T^ . Specifically, 
dN3/dy = dA^'=''+"™*7dy = 2800. Assuming local ther- 
mal equilibrium, for a gluon-dominated plasma we find: 



T(t,x)= V^V(t,x)/16C(3) ,r >ro 



(13) 



The Debye screening scale is then given by = AiragT'^ 
and the relevant gluon mean free path is easily evaluated: 
Xg = 1/crSffp with = {9/2)T:al/ml,. 

In our simulation we generate in-plane jets, 0jet — 
^reaction piano = 0, sincc in Central Pb-I-Pb collisions 
the interaction region is nearly azimuthally symmetric. 
We evaluate the differential bremsstrahlung spectrum 
Eq. (jlOp to the lowest non-trivial interference between 
the gluon emission from the hard jet production and the 
subsequent interactions in the QGP [l9|. In our simula- 
tion the allowed gluon phase space is Aqcd < uj < E, 
Aqcd < k± < 2uj [40| . and the transverse momentum 
transfers between the jet and the medium are in the in- 
terval < q±i < \fsji = ^JmDE/2 [li]. Finally, we 
note that we study two coupling strengths between the 
fast partons and the medium, = 0.3, 0.5, and use a 
running as{kT) for the emission vertex. 



IV. PREDICTIONS FOR HEAVY-ION 
COLLISIONS AT LHC ENERGIES 

With the numerically expensive simulations of en- 
ergetic quark and gluon propagation through the ex- 
panding QGP medium created in central Pb+Pb reac- 
tions at the LHC completed, we proceed to evaluate the 
Z°/7*-|-jet cross section in these ultrarelativistic nuclear 
collisions. In our approach, partonic energy loss is in- 
corporated through two crucial ingredients. The first is 
the fraction of lost energy that falls within a certain cone 



radius R and is carried away by gluons of energy above a 
minimum energy cut-off LUram |10l |. We denote this frac- 
tion by: 



E) 



AE{R,Lo^in,E) 
AE{R°°,0,E) 



(14) 



In Eq. E is the initial energy of the jet parent par- 
ton, and the subscript q, g indicates whether this parton 
is a quark or a gluon. The function fq^g describes the 
fraction of partonic energy loss that will be redistributed 
inside the jet for fixed jet reconstruction parameters R 
and Wmin- A_B(i?, cjinin, _E) is obtained from the double 
differential spectrum of medium-induced bremsstrahlung 
dl^ /dujdr as follows: 



AE{R,uj,r,in,E) = / duj dr 

Juitnia Jo 



dis 
diddr 



{u,r). (15) 



Figure [S] shows an example of the mean 
AE{R,u;min,E)/E for 100 GeV quark and gluon 
jets that were created in central Pb-f-Pb reactions at the 
LHC and have propagated through the QGP medium. 

The second ingredient we will need is Pq.g{e), which is 
the probability distribution for the fractional energy loss 
of an energetic parton due to multiple gluon emission, e ~ 
'Y^^ujijE. The evaluation of this distribution from the 
calculated medium-induced bremsstrahlung spectrum is 
described elsewhere [s^]. Knowledge of Pq,g{e) is essen- 
tial for the probabilistic superposition of vacuum and 
medium-induced parton showers and the calculation of 
inclusive and tagged jet cross sections in heavy in reac- 
tions. 
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FIG. 5: (Color online) Fractional energy /S.E/E lost by 
100 GeV quark and gluon jets inside a cone of radius R via 
gluon emission of a; > cJmin as a result of final-state interac- 
tions in the QGP created in central Pb-|-Pb reactions at the 
LHG. 



A. Evaluating the suppressed ^"/7*+jet cross 
section 

Consider the invariant cross section for Z/7*-tagged 
jets in p-l-p collisions at an arbitrary order in perturba- 
tion theory, which we write as: 



da 



!/(z) 

y(z)„ 



dy(2 



y(jot) , 

!/(jot) m 

da 



(16) 



jot) 



Here, we are interested in obtaining the medium-modified 
jet cross section per binary N + N scattering in terms of 
the above expression. In particular, medium-induced en- 
ergy loss docs not affect on average the direction of the 
jet (and y(jot)) a-nd its effects arc contained primarily in 
P(jet)- We define pg P(jot)(l - (1 - fq.g)<^) which is 
the reduced momentum accounting for medium-induced 
energy loss that falls outside a certain cone radius and 
energy cut-off. In the limit of zero energy loss (e — > 0) 
or very large jet cone radius (/g,g — > 1) one recovers 
Pq — P(jct)- Accounting for the fact that quark and 
gluon fractional energy loss is described by a probability 
distribution Pq,g{() and the transverse momentum Jaco- 
bian: 

d^V>{iet) 



d2p 



(Q) 



[l-(l-/,,,(a;„i„,i?))e]2 ' 



for the quenched jet cross section we find: 
da ^ r , „ , . 1 



= E/ deP.A^) 



X ■ 



[I - {I - fg,giUrain, 



da'i^'J 



PQ 



rf^P(Z)rf2p(jct) V[l - (1 - /g.ff)^)] 



The physical meaning of Eq. pT)) is that the observed 
tagged jet cross section in A-l-A reactions is a proba- 
bilistic superposition of cross sections for jets of higher 
initial transverse energy. This excess energy is then 
redistributed outside of the jet due to strong final- 
state interactions. Here, Pq,g{e) and fq,g(ujmin, R) are 
obtained using the GLV formalism for evaluating the 
medium- induced gluon brcmsstrahlung [l9l [25j . Also, 
da'^'^ /d'^P(^z)d'^P(jct) are the differential cross sections for 
away-side quark and gluon jets, respectively. We calcu- 
late the relative fraction of these jets to lowest order in 
perturbation theory. Even though at next-to-leading or- 
der one might expect a small correction to these results, 
we show in Appendix [Q that a variation of the 0{as) m 
the relative quark to gluon jet fraction does not affect 
the predicted quenching of Z°/7*-tagged jets. 

It is often the case that the Z°/7*+jet final-state chan- 
nel is measured without placing restrictions on the mo- 
mentum of the vector boson 0]. In this case, integrating 
over P{z} in the LHS and RHS of Eq. ([TT]) we obtain: 



(17) 



da 
rf^PQ 



1 



[1 - (1 - /g,3(Wmin,^))e 



Pq 



da'i'S 



d^P(jct) \[1 - (1 - Ug)e)] 



(18) 



For all practical purposes, Eq. ([T8|) contains the same 
physics as the suppression of the cross section in inclusive 
jet measurements [l^, [13, HO] ■ We will demonstrate this 
shortly in our numerical results section. 

In the special case of a tree level calculation for the 
elementary quark and gluon jet tagged cross sections we 
recognize that: 



da 



da 



H J2 ~ J2 ^ (P(^) ^ P(jct)) ■ 

Substituting this result in Eq. (|17p , after straightforward 
algebraic manipulation wc find: 



da 



^{4>(z) " 0(jct) - tt) da'i-3 
d^V[z)d^VQ ^ Vt {q)Vt {Z){^ ~ fq,g) d^P(z) 



E 



1 - Pt iQ) / Pt (Z) 
1 - Us 



(19) 



Note that in Eq. (jl9p the remaining (5-function sim- 
ply reflects the fact that at tree level the Z°/7* and 
the jet are exactly back-to-back and the angular dis- 
tributions can be trivially integrated over. From the 
< fq.gi'^mini R) < 1 and the properties of Pq,g{e) one 
finds that the double differential cross section exists in 
the region fq.gPr (z) < Pt (q) < Pt {z) ■ 

It is instructive to investigate the fq^g — > 1 limit in 
Eq. (Ull). For pr (q) < Pt{z), PqA^ > 1) = en- 
sures a lack of singularity. If Pt(Q) ~^ Pt(z) we have 
lim,^oA,s(e) ^ exp(-(iVf),,g)(S(e) [H. Here, (7VS),.g 
is the average number of medium-induced gluons for a 
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quark or gluon jet respectively. To interpret the result- 
ing cross section we integrate over (q) to find: 



da 



d2p 



Pt (Q) =Pt (z) 



-iN^K.^ (20) 



d2p 



The physical meaning of Eq. ([20)) is that at lowest or- 
der the coincident Z'^/j*+jet production with Pt(q) = 
— Pt is directly proportional to the probability not to 
lose energy via medium-induced gluon bremsstrahlung. 
Eqs. ([T5)) and (PO)) imply that one might gain easy access 
to {N^)q^g and Pq^g{e). However, higher order correc- 
tions significantly alter these relations. To disentangle 
these quantities/distributions from tagged jet measure- 
ments one still needs accurate baseline calculations of the 
selected final-state channel at NLO and, ideally, precise 
experimental measurements. 

We conclude this section by restating in a compact 
form the basic results for the evaluation of inclusive and 
tagged jets in the ambiance of strongly-interacting mat- 
ter: 



da 



dp' 



'T{Q) 



dpT 



(jot) 



(21) 



da 



dpT(z)dpT 



{Z)"-PT(Q) 



dai'9 



dpT {Z)dpT iiet) 



(22) 



B. Numerical Results 

We now present numerical results for the QGP- 
modified Z'^/7*-|-jet production cross sections in cen- 
tral Pb-hPb collisions at the LHC at ^/snn = 4 TeV. 
As for p+p reactions, we consider jets tagged with 
Z/j* fi^ + fi~ using acceptance cuts of \y\ < 2.5 
for both jets and final-state muons. Also, the invariant 
mass of the muon pair is required to be in the interval 
mz ± 3r^, where mz = 91.2 GeV and = 2.5 GeV. 
In this paper we focus on the effect of final-state interac- 
tions, motivated by the fact that both the large virtuality 
= ?^^T(z) ^-'^'^ range of mz < rriT (z) ^ \/snn/2 
(at midrapidity) are expected to reduce the significance 
of initial-state effects relative to RHIC kinematics [l3 |. 
We defer detailed studies of initial-state interactions to 
future work. 

We start with the tree level results for the px- 
differential jet cross section where the tagging Z/j* is 
required to have 7 GeV < pt < 13 GeV and 92.5 GeV < 
Pt < 112.5 GeV. These cross sections are shown in the 
top and bottom panels of Figure [51 respectively. In the 
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FIG. 6: (Color online) LO results for the-pr differential 
Z/7*-|-jet cross section where the tagging Zj^* is required 
to have 7 GeV < pr < 13 GeV (top panel) and 92.5 GeV < 
■pT < 112.5 GeV (bottom panel). We present simulations for 
y/SNN = 4 TeV p-|-p and central Pb-|-Pb collisions at the LHC 
for three different jet cone radii R. 



collinear factorization approach at leading order the mo- 
mentum of the Z^ /j* is exactly balanced by the mo- 
mentum of the jet. Furthermore, to this order the jet is 
represented by the energetic parent parton and the cross 
section is insensitive to the choice of the jet cone radius 
R [4l| . In Pb-|-Pb reactions the redistribution of the jet 
energy is determined by the medium modification to the 
parton shower but the experimentally observable effect 
is controlled by the jet cone radius R. At tree level, the 
choice of radius (and/or Wmin) also determines the px 
range where the tagged jet cross section does not vanish, 
see Eq. We note that, unless stated otherwise, in 

our numerical simulations we assume no momentum cut 
on the soft particles/partons, Wmin = GeV. The dis- 
tribution of jets in central Pb-|-Pb collisions in Figure [6] 
reflects the probability distribution for the fractional en- 
ergy loss of partons (^'(e)) [3l|; the average being over 
the quark and gluon jets. 
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The inserts in Figure [6] show a measure of nuclear ef- 
fects for tagged jet cross sections in heavy ion reactions 
relative to the binary collision (A'^bin) scaled cross section 
in proton-proton collisions: 



10^ 



1 



da 



AA 



dan 



(iVbin) dpT{Z)dpT(Q) 



dpT(Z)dpT (jot) 

(23) 

As we will see shortly, is much more useful at NLO. 
At tree level it is limited to the /^* trigger momen- 
tum range. The large suppression toward the upper edge 
of transverse momentum interval is reminiscent of the 
exp(— (iVS)g_g) factor in Eq. (|20|) and indicates that par- 
ton energy loss in the QGP proceeds through the emis- 
sion of several semi-hard (of energies of the order of 
few GeV) gluons. Unfortunately, in realistic measure- 
ments the finiteness of the pt trigger bin itself will pre- 
vent a more accurate determination of {N^)q,g. 

The NLO pT-diffcrential jet cross section with a tag- 
ging Z/^* required to have 7 GeV < < 13 GeV 
is shown in Figure [T] This is the first order in per- 
turbation theory at which the distribution of jets rela- 
tive to the constrained Z /^* momentum can be studied. 
We present simulations for three different values of the 
jet cone radius, R = 0.2, 0.4, 0.8, and two different val- 
ues for the jet-medium coupling strength, = 0.3,0.5. 
The main physics effect illustrated by this figure is the 
QGP-induccd modification to the vacuum parton shower. 
Specifically, its broadening HEj implies that part of 
the jet energy is redistributed outside of the jet cone 
and the differential jet distribution is downshifted toward 
smaller transverse momenta. The smaller the jet cone 
radius the more pronounced this effect is, as is clearly 
seen in the panels of Figure [71 As the jet cone radius 
is increased, the medium-modified curves approach the 
p-|-p result, as more and more of the medium- induced 
bremsstrahlung is recovered in the jet. The variation in 
the magnitude of at transverse momenta larger than 
the Pt of the tag is controlled by the shape of the jet spec- 
trum. The most striking feature is the sharp transition 
from tagged jet suppression above pt (z) to tagged jet 
enhancement below Pt (z)- This transition is a unique 
prediction of jet quenching for tagged jets and will be 
experimental evidence for strong final-state interactions 
and parton energy loss in the QGP. Note that while the 
same energy redistribution occurs for inclusive jets , 
the monotonically falling spectrum prevents the observa- 
tion of such enhancement in inclusive jet measurements. 

Next, we consider the same cross sections and ratios 
for the tagging Z^ /^* now required to have 92.5 GeV < 
Pt < 112.5 GeV. Our results are shown in Figure[8] The 
cross section enhancement in Pb-|-Pb relative to p+p for 
Pt {jet) ^ 100 GeV is quite dramatic. This is particularly 
true for the smallest cone radius, R = 0.2, where the 
medium modified result may become an order of magni- 
tude larger than in p-l-p, see Iaa in the inset of the upper 
panel. The fiuctuations in P^a ^-rise from the numerical 
accuracy of the simulation and the rigid acceptance cuts 
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FIG. 7: (Color online) The NLO pr-differential cross section 
per nucleon pair for jets tagged with Z /'y* — >■ jj.'^ + fj,~ in p-|-p 
and central Pb+Pb reactions. The tagging Z/^* is required 
to have 7 GeV < pt < 13 GeV. Results are shown for three 
different values of jet cone radius R — 0.2 (upper panel) R = 
0.4 (middle panel) and R = 0.8 (lower panel). The ratio of 
the tagged QGP-modified cross section per nucleon pair to 
that in p+p, Iaa, is shown in the inset of each plot. 
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FIG. 8: (Color online) The same cross sections and ratios are 
shown as in Figure [7] but now the tagging Z/y* is required 
to have 92.5 GeV < pr < 112.5 GeV. Notice the sensitivity 
to the jet cone radius as R varies from 0.2 to 0.8 (upper to 
lower panels) and also the significant enhancement in Pb+Pb 
relative to p+p for PT{jct) lS 100 GeV. 



for both decay muons. In practice, these fluctuations "will 
be smeared out. We conclude that a transverse momen- 
tum for the Z'^/j* tag on the order of pT ^ nT-z will 
be optimal to establish the QGP effects on jet propa- 
gation and distortion in heavy ion reactions if sufficient 
integrated luminosity becomes available. 

We finally present results for jets tagged -with Z°/7* — >■ 
/i"*" -|- /Li~ without imposing restrictions on the transverse 
momentum of the tagging particle. This is precisely 
the way in which the Tevatron II measurements were 
done Q . The resulting cross sections are essentially cross 
sections for a special subset of processes that contribute 
to inclusive jet production. The alteration of these "in- 
clusive" cross sections in high-energy nucleus-nucleus col- 
lisions can be studied through the nuclear modification 
ratio: 



R'^AAi-^^ '^min) — 



dcTAA 



da, 



pp 



(Nun) dpT (Q) / dpT (jet) 



(24) 



The NLO py-difFercntial result for Raa is shown in Fig- 
ure [5] In the top panel we fixed Wmin = GeV and varied 
R. Bands represent a variation in the coupling strength 
of the jet parent parton to the QGP medium between 
as = 0.3 and as = 0.5. The strong dependence on the jet 
radius R emphasizes the pronounced difference in the an- 
gular distribution of vacuum and medium-induced parton 
showers in the framework of the GLV approach to par- 
ton energy loss in the QGP [35|. In the bottom panel of 
Figure [S] we fixed R — OA and varied Wmin- The observed 
additional suppression in Ra\ reflects the additional en- 
ergy lost through the elimination of the medium-induced 
semi-hard partons of px < Wmin from the jet definition. 
The insert shows that there is little difference in the sim- 
ulated Raa '^min = GeV and Wmin = 2 GeV and, 
consequently, that the medium-induced gluon spectrum 
in the heavy ion reactions at the LHC are typically of 
larger energy (oj ^ 5 GeV). 



V. CONCLUSIONS 

Pb+Pb collisions at the LHC at ccnter-of-mass ener- 
gies up to 5.5 TeV per nucleon pair are expected to con- 
solidate the evidence for the creation of the QGP in ul- 
trarelativistic heavy ion reactions. Such collisions will 
also probe the plasma properties in a new temperature 
and energy density regimes and open final-state channels 
for jet tomography [ij of the QGP that are currently 
inaccessible to experiment or limited by statistics. One 
of these channels, and the focus of the work presented 
in this paper, is the Z'^/7*+jet measured through the 
Z'^/l* — >■ M"'' + dilepton decay channel. Because of 
the short Z° production time and the fact that the decay 
dileptons reach the detectors unscathed by the strong in- 
teractions, this channel has been proposed as a way to ex- 
perimentally constrain precisely the initial energy of the 
associated jet and through comparison to theory quantify 
the jet quenching properties of the QGP |ll| - [l3| . 
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FIG. 9: (Color online) NLO results for jets tagged with 
Z/7* — >■ presented as the ratio of the QGP-modified 

cross section per nucleon pair to that in p+p. Note that 
there is no pt restriction on the Zj-^* . In the top panel we 
fix Wmin = and vary R, whereas in the bottom panel we fix 
R = 0.4 and vary cjmin. The insert shows that the variation 
of 7?AA(<^min) for Wmin < 5 GeV is negligible. 



In this paper we examined the tagging power of the 
Z"/7* in p+p and A+A collisions at leading and next- 
to-leading orders. The main advantage of the NLO cal- 
culation is the ability to precisely predict the transverse 
momentum distribution of jets associated with a dimuon 
tag in a narrow interval. We found that, in contrast 
to the naive tree-level result, to O{Gf0?^ there is very 
large ^ ±25% uncertainty in the jet transverse energy 
distribution relative to the dimuon tag pT- Our results 
demonstrate that direct experimental determination of 
the jet fragment distribution Pt (hj/PT (jet) without ade- 
quate NLO theory is not realistic. Some of the variation 
in the transverse momentum of the jet is likely to be 
eliminated by additional kinematic cuts but at the cost 
of reducing the already small cross section. We conclude 
that it is not clear whether jet tagging with Z°/7* at the 
expense of a substantially reduced production rate offers 



significant advantages relative to direct photon tagging. 

We also presented for the first time a calculation of 
tagged jet cross sections in heavy ion collisions up to 
next-to-leading order. Together with earlier tree-level 
and one-loop results for inclusive jets [l^ [l^ , our work 
paves the way for more detailed and systematic theoret- 
ical and experimental studies of jet propagation in the 
background of strongly-interacting non-Abelian plasmas. 
It builds up a suite of differential jet observables to probe 
the microphysics behind the in-medium modification of 
parton showers. This new direction of research is promis- 
ing in that it has the potential to overcome the limitations 
of leading particles [22| and leading particle correlations 
as tomographic probes of the QGP. Comparison between 
jet theory and upcoming experimental measurements at 
RHIC and at the LHC will help constrain theoretical 
model assumptions and approximations to parton prop- 
agation in hot and dense nuclear matter. Elimination of 
unreliable phenomenology will, in turn, dramatically re- 
duce the systematic error bars on the extracted plasma 
properties. 

As a first example of tagged jet cross sections in nuclear 
collisions we calculated Z''/7*-fjet production process in 
^/sj^ = 4 TeV central Pb-t-Pb reactions at the LHC. We 
demonstrated that without constraints on the momen- 
tum of the Z'-' /^y* decay dimuons, the quenching of tagged 
jets as a function of the jet cone radius R and momentum 
acceptance cut Wmin exhibits all features characteristics 
of inclusive jet measurements. Specifically, the variable 
attenuation rate of jet production R''^\{R, uj-min) provides 
first-hand information about the differential spectrum of 
the QGP-induced gluon bremsstrahlung. In contrast, 
/^"^^(i?, Wniin) yields additional insights into the charac- 
teristics of the medium- induced parton showers. We have 
found that direct determination of the mean number of 
medium-induced gluons [sif is hampered by NLO correc- 
tions. However, the dramatic transition from suppression 
of tagged jets above the Z'^ /j* px to strong enhancement 
below this reference transverse momentum characterizes 
the probabilistic distribution of the energy lost by ener- 
getic partons as they traverse the QGP medium. Exper- 
imental measurements of such jet enhancement will pro- 
vide an unambiguous confirmation of the dominant role 
of final-state interactions in the experimentally observed 
modification of jet and particle production cross sections 
in nucleus-nucleus collisions relative to an independent 
superposition of nuclcon-nucleon scatterings. 

We conclude by emphasizing that one can also take 
advantage of tagged jet physics in heavy ion collisions 
by using alternative final states. In particular, direct 
photon-tagged jets are characterized by significantly 
larger cross sections when compared to Z''/7*-tagged 
jets and are accessible in central Au-|-Au collisions at 
RHIC. It is a top priority to provide accurate theoretical 
predictions [ij] for these reactions and to study how the 
QGP responds to propagating hard partons [3^ tagged 
by an electrowcak boson. 
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Appendix A: Lowest Order cross section for 
2'''/7*-tagged jets 

The lowest order cross section for the ^°/7*-l-jet pro- 
duction in hadronic collisions is well-known Never- 



theless, it is useful to outline a few selected steps in its 
derivation in order to clarify the kinematics and produc- 
tion rate for this process. To LO, only four diagrams, 
given in Figure I10[ contribute. Our convention is to la- 
bel the momenta as given there. Specifically, we have 



pi = (a;iPi,0, 0,xiPi) , 

P2 = (2^2^2,0, 0,-0:2^2) , 

P3 = (mT cosh 2/3, pT,mTsinh 2/3 ) , 

P4 = (j>T cosh 2/4 , - pt , Pt sinh 2/4 ) , 



where in all cases P3 refers to the and P4 to the final- 
state quark or gluon jet. Additionally, Pi = -P2 = 
xi , X2 are the momentum fractions carried by the initial- 
state partons, and rriT ~ ^p|, + m^. 

The amplitude for gluon jet production (upper panel of 
Figure [To)) is the sum of two matrix elements Mgi +Mg2. 
Specifically, for the q + q ^ Z° + g process we have: 



iMgi 

lMg2 



■ 9sg 

I — 



.(P2, S2){YTy-i{r) {R, (1 + 7') + Lg{l - 7')) SlXe* , 

V2 Pa 



■ 9sgz 



3(P2, S2)(7^) {R, (1 + 7') + Lg{\ - 7')) ^{YT^)u{pus^)elel 



Pb 



(Al) 
(A2) 



where pa = P2~Pi and p^ = Pi — P4. The contribution for quark jet production (lower panel of Figure [TU]) is the sum 
of Mqi + Mq2 for the process q + g —i' Z'^ + q: 



■ 9s9z 
* — ^ 
V2 



(P4, s,){YT'^)Jf{r) {R, (1 + 7') + L,{1 7')) u{p,,s,)e^el , 
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Pa 



■T^f -989. 

v2 



(P4, S4)(7^) {R, (1 + 7') + - 7')) %(7''r'^)"(Pl, Sl)6.6; 

Vc 



(A3) 
(A4) 



r 



where Pc = Pi +P2- In what follows, we use the notation 
of Ref. [131 and have defined: 

(A5) 



23^^(1 ^-u;) 

with Xw = sin^ « 0.223 for the on-shell Weinberg 
angle and a « 1/137.04 [2^. Also, we have introduced: 

i?2 = 4e2 sin4 0,„ , 

L\ ^tI~\ eg Tq sin^ 6*^ + 4 sin" 6*^ , 

where is the weak isospin of quark q (that is, r = 1 
for u, c, t and r = — 1 for d, s, b) and Cq is the fractional 
electric charge of quark q. 



In obtaining the squared amplitudes, |Afgp = \Mgi 4- 
A/g2p and \Mq\'^ = \Mqi + A/g2p, we perform the ex- 
ternal gluon polarization sums with the replacement 
^tiv — >■ —9fj.u- For massive vector bosons these sums 
read: 



9p.u 



{p-i)ti{p?.)i' 



(A6) 



After straightforward color and Dirac algebra manipula- 
tions the final tree-level result is obtained as: 



\M„ 



st 



, (A7) 
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FIG. 10: LO diagrams contributing to the Z^+jet production. 
Here, p3 is the momentum of the Z". At LO the contribution 
from quark and gluon jets are conveniently separable; how- 
ever, at NLO real gluon emission complicates that distinction. 



where Us ~ <7s/47r and the Lorcntz- invariant Mandel- 
staam variables are: 



where xi and X2 are the solutions to the coupled equa- 
tions: 

(xi -t- a;2)\/5/2 — cosh 2/3 — pt'4 cosh 11/4 = 0, 

(xi — X2)^/S/2 — sinhy3 — pj-^ sinhy4 = . 
We find: 



Xi 



X2 



Vs 



(All) 
(A12) 



and LO cross section is then given by (j)t ~ Pt^ ~ Pt^) 



"^"l^ = E J^}\ ■ (A13) 



dys dy. 



s = {pi +P2f , i= [pi -ps) 



{P2 -Paf 



The 7*-tagged jet matrix elements squared can be ob- 
tained from Eqs. (|A7I) and (IA8|) with the replacement: 



(A9) 



as well as letting rn^ — > Q^, where Q is the photon vir- 
tuality. While in p-t-p reactions we are interested in the 
total jet rate production rate, oc |Mgp + |Mgp + |M^p 



the separation into final-state quark 



and gluon jets will prove useful in the treatment of the 
medium-induced radiative energy loss in A-l-A reactions. 
Note that when the initial-state partons are interchanged 
the matrix elements can be trivially obtained with the 
substitution i -H> u. 

In the collinear factorization approach the differential 
■^°/7*+jct cross section is evaluated as follows: 



da 



E 



dxidx2 f{xi,Q)f{x2,Q) 



X {2nf6Hpi 



■P2 -Ps-Pi) 



2X1X2S 

d^P3 d^p4 
2£;3(27r)3 2£;4(27r)3 ' 

(AlO) 



Further, one notices that d^pi ~ Ei dyi d'^pT^ and 



S'^{pi+P2-P3-P4:) = ■^S{xi-xi)6(x2~X2)S'^{pts-PtJ , 



Appendix B: The Z/j* 1+ + T Dalitz Decay 

In this Appendix we discuss the /j* decay to dilep- 
tons with the goal of addressing the experimental mea- 
surements at the Tevatron and at the LHC. Our approach 
is to work in the rest frame of the parent particle. If this 
particle has mass mz-, its differential decay rate is given 
by: 



dT 

an 



327r2 



\M 



,2 IPI 



(Bl) 



where p is the momentum of one of the final-state parti- 
cles (the other particle has momentum — p) and |Mp is 
the matrix element squared for this decay. For instance, 
in the case of 7* p^ + fi^: 



iM 



-'ieu{p,S2)Yv{p,si)el 



3^ 

sp'in 



\M\ 



(B2) 



where Q is the virtuality of the photon (that is, the in- 
variant mass of the muon pair). The calculation of |A/p 
for Z ^ + p,~ is analogous: 



iM = -i^^iP, ^2)7" (^e(l + 7') + L,il - 7')) v{P, si)e: 



^Ei^i^ 



4glml{Rl+Ll) 



(B3) 



spin 



14 



where Re = 2xw and = Ixyj — 1 (sec Appendix El for 
the definition of and gz ) • Clearly, the distribution of 
decay products is isotropic. 

To evaluate the cross section for the Z°/7*(— >■ + 
/^~)+jet process we use the ex per imentally determined 
branching ratio F^^/F w 0.037 [2g. The velocity of the 
boson in the lab frame, in terms of the variables dis- 
cussed in Appendix [Aj is given by: 



/3 



(pT,mT sinh ^3) 
m,T cosh 7/3 



(B4) 



In the Z'' rest frame, the 4 momenta of the two decay 
muons are given by: 



Viz 



mz 



(1, sin 6' cos 0, sin sin (f), cos 6) 



^(1 
2 ^ ' 



where 6 and </> are uniformly distributed over the solid 
angle. The 4 momenta of the two decay muons in the 
lab frame can be obtained with an appropriate Lorentz 
boost PL = Upz where U is given by: 



/ 7 














\ 


/3.7 


1 + 


(7-l)ff 


(7 




(7 


- 




Pvl 


(7 




1 + 


(7-1)1 


(7 


^)-^ 




V/3^7 


(7 


- 1)^ 


(7 




1 + 


(7-l)fl 


/ 



This approach enables one to determine the distribution 
of the final-state dileptons and to simulate any relevant 
experimental acceptance cuts. 



Appendix C: The role of quark versus gluon energy 
loss in inclusive and tagged jet quenching 
measurements 

At next-to-leading order perturbative splitting pro- 
cesses [q{q) q{q)g, q{q) gq{q), g gg, g qq] 
complicate the flavor identification of the jet parent par- 
ton. In proton-proton collisions the ambiguity of quark 
vs gluon jets primarily affects fiavor-tagged jet observ- 
ables d^. In nucleus- nucleus collisions the primary con- 
cern is the possible uncertainty in the amount of medium- 
induced gluon bremsstrahlung, which is proportional to 
the average squared color charge {Cp =4/3 and Ca = 3 
for quarks and gluons, respectively). When parton en- 
ergy loss is not very small relative to the initial jet energy, 
the factor of Ca/Cf = 2.25 difference between Ai?g and 
Ai?g is greatly reduced (s^. It is nevertheless instruc- 
tive to examine how changes in the relative abundances 
of quark and gluon jets (possibly arising in higher order 
calculations) affect the observable R^aa ^'^'^ ^^aa- 

A solution that we use for inclusive [l3| and tagged 
jets is to evaluate at tree level the quark and gluon jet 
fractions Uq, Ug [uq + iig = 1) for the specific final state. 
In this appendix we demonstrate that the sensitivity of 



T3 



D 
T3 
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+ n /n„(1 +/ 


aJ2) variation 


. _ Gluon fraction 


+ n /n„(1 +/- 


aJ2) variation 


. "q + "g = 1 s^'^ = 4 TeV 















100 



p/*[GeV] 



200 



300 



FIG. 11: Quark and gluon jet fractions Uq, Ug versus pr 
for Z''/7*-tagged jets in ^snn = 4 TeV p+p collisions at the 
LHC (evaluated at tree level). The band represents a possible 
variation of Ug/uq of ±as(mT z) /'i- 



inclusive and tagged jet suppression measures Raa and 
Iaa to 0{as) corrections to this quark versus gluon jet 
mixing is small. Parton splitting is proportional to the 
strong coupling constant and we focus on a possible vari- 
ation in Ugjuq by this amount. Figure [TT] shows the 
calculated quark and gluon jet fractions: 



dp7 



dp7 



d-pi 



dpi 



(CI) 



in y/sNN ~ 4 TeV p-f-p collisions at the LHC. The de- 
tailed behavior of these fractions as a function of px, 
including the non-monotonic structure around pt ~ mz, 
arises from the interplay of the initial-state abundance of 
(anti)quarks and gluons and the behaviour of the short- 
distance scattering, Eqs. (|A7p and (|A8I) . Bands represent 
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FIG. 12: Effect of Uq, Ug variation on the suppression of the 
Z''/7*(— >■ /.i^ -|- /x~ )-tagged jet cross section in central Pb4-Pb 
reactions at ^/sNN = 4 TeV at the LHC. We have chosen a 
jet cone radius _R = 0.2. 
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Due to the smallness of strong coupling constant at large 
Pt, the size of the bands in Figure [TT] is also small. 



The effect of such Uq/ug uncertainty on is shown in 
Figure [12] for a dimuon tag 92.5 GeV < pt < 112.5 GeV 
and a jet cone radius R = 0.2. Unlike in Figure [3 we 
have used a much smaller transverse momentum range 
for the jet and a linear scale to illustrate that the 
resulting uncertainty in the quenching of Z'^/7*-taggcd 
jets is of the order of 1%. 



FIG. 13: Effect of Uq, rig variation on the suppression of 
the inclusive jet cross section in central Au+Au reactions at 
y/SNN = 200 GeV at RHIC. We have chosen a jet radius 
7? — 0.4. Insert shows the same effect for the ratio of two 
cross sections with R = 0.2 and R = 0.4, respectively. 



the following variation of the gluon to quark jet ratio: 



1± 



asirriTz) 



(C2) 



Finally, we point out that for inclusive jets the vari- 
ation in the cross section suppression i?^^^ is similarly 
small. The is illustrated in Figure [T3] for RHIC ener- 
gies of y^SNN = 200 GeV. Our calculation is for central 
Au-|-Au reactions with a choice for the jet cone radius 
R = 0.4 [131 • Everywhere in the accessible pt range the 
uncertainty on the attenuation of the jet rate in the final 
state is less than 5%. 
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